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Abstract 

The remagnetization dynamics of monolayer dot array superlattice XY 2-D spin model with 
dipole-dipole interactions is simulated. Within the proposed model of array, the square dots are 
described by the spatially modulated exchange-couplings. The dipole-dipole interactions are ap- 
proximated by the hierarchical sums and spin dynamics is considered in regime of the Landau- 
Lifshitz equation. The simulation of reversal for 40000 spins exhibits formation of nonuniform 
intra-dot configurations with nonlinear wave/ anti- wave pairs developed at intra-dot and inter-dot 
scales. Several geometric and parametric dependences are calculated and compared with oversim- 
plified four-spin model of reversal. The role of initial conditions and the occurrence of coherent 
rotation mode is also investigated. The emphasis is on the classification of intra-dot or inter- 
dot (interfacial) magnetic configurations done by adaptive neural network with varying number of 
neurons. 

PACS: Nanostructures: 79. 60. J, Magnetostatics 41.20.G, Magnetic ordering 75.10, Numerical simula- 
tion studies 75.40.Mg 
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1 Introduction 

Many interesting applications of the nanometer periodic magnetic dot arrays are expected in the 
magnetic recording [1] and magnetic sensors. The magnetic nanostructures represent model systems 
for study of interesting physical phenomena stemming from the well controlled dot shape anisotropy 
and dominance of magnetostatic interactions at the scales comparable or larger than dot size. The 
magnetostatics cause the coupling of domains over the inter-dot interfaces. The current perspective 
technologies require understanding and controlling of dynamic magnetization processes in fine particles 
on nanoscales, because by them the magnetic bits are written into recording media. Hence of particular 
interest is the dynamics of magnetization reversal. 

The static and dynamic nonuniform magnetization states in magnetic nanostructures including 
the thin films and isolated small particles have been studied by micromagnetic [2, 3, 4] and Monte 
Carlo [5, 6, 7, 8] simulations. On a qualitative level intra-dot nonuniformities can be characterized as 
a superpositions of vortices, fiowers, domain walls, etc.. The problem of nonuniformities in exchange- 
dominated structures can be solved exactly [9]. 

Our recent activity is directed to algorithms allowing systematic and fast classification of simulated 
magnetic structures. The methodology presented in this paper is inspired by the progress in the 
theory of artificial neural networks that are nonlinear models suitable to reduce, classify or interpolate 
data structures representing the complex patterns [10]. In the previous paper [11] we suggested the 
implementation of radial basis function networks to model the ordering of dots in arrays. In the next 
paper [12] the neural network training was applied to multiscale numerical computations. The lack 
of our previous two approaches was the absence of realistic dynamics and restriction to quasi-static 
limit. In the present paper our interest is focused to the networks which analyze spin configurations 



generated by the spin lattice dynamical simulations. It means that network solves here classical post- 
processing problem known in the micromagnetic context [13]. Our approach has many in common 
with experiment treatment [14], where self-organized network was employed to identification of vortex 
modes of magnetic thin film media. In the presented paper the classification is applied to configurations 
generated by the magnetization reversal within the XY superlattice model of ultra-dense array. The 
motivation for study of ultra-dense spin structures was the expectation of extremal inter-dot effects. 

The main reason for the choice of the XY spin model is relative computational simplicity of 2-D 
dynamics. But there are two regimes for which 2-D approximation is justified physically, (i) The first 
regime belongs to limit of very strong damping which suppresses precession and accelerates reversal 
[15]. (ii) The second regime represents limit of strong in-plain [16] or shape anisotropy associated with 
demagnetization field of thin film sample [17]. In this regime anisotropy causes deformation of spin 
trajectory and accelerates reversal [18]. Nevertheless, one should to note that realistic damping is 
sufficiently small, and thus 2-D remagnetization is practically quasistatic. In the quasistatic limit the 
damping parameter has meaning of auxiliary parameter rescaling the time. Recent results obtained for 
discrete spin systems provide more optimistic view on the valuation of artificial dynamics of reversal. 
The example is the simulation of [19] demonstrating that switching time Monte Carlo dependences 
can coincide after some rescaling with predictions of the stochastic Landau-Lifshitz equation. 

The plan of the paper is the following. The interaction picture and dynamics of array model are 
defined in section 2. The simulations of array reversal are presented and discussed in section 3. The 
details of adaptive neural network treatment of intra-dot configurations are given in subsection 3.3. 

2 XY periodically modulated superlattice model of array 

In this section we formulate the interaction picture and dynamics of ultra dense array model. The 
magnetic ordering is described by the systems of the classical planar spins located on the square two 
dimensional lattice. The spin on site a is affected by effective field 



where Sq, = Sax^x + Say^y is the planar spin, H is the external magnetic field in energy units, L is the 
linear size of the square L x L spin lattice system. Here and in further Cx, Gy denote the unit vectors 
along the main Cartesian axes. Rather convenient for a computer purposes is to use one dimensional 
enumeration scheme of spin sites a G {0,1,2,...L^ — 1}. (All of the geometric details presented in this 
section are discussed in the Appendix A.l.) The first sum of Eq.(l) represents dipolar contribution 
proportional to constant D, where r^,^ is vector pointing from site a to (3. The vector is measured in 
lattice constant units, n^^ = Vap/lrapl is the auxiliary unit vector. The second sum of Eq.(l) involves 
nearest neighbor exchange coupling term Jeag determined by single parameter J > multiplying the 
site-site connectivity matrix e^^ G {1)0} constrained by \ra0\ = 1 condition; e^^ is for any bond at 
inter-dot boundary and 1 inside a dot. The site dependence of e^/j takes into account square shape 
of dot of the linear size Lq. The spin system partitioned by e^/j consists of = (L/Lq)^ dots with 
minimum inter-dot distance equal to the unit lattice spacing and corresponding to extremal space 
filling. 

In 2-D case Landau-Lifshitz equation reduces to damping term. The numerical integration scheme 
discussed in Appendix A. 2 can be derived from recursive propagator form 



where ^^(r) = (Sa(r) • ex)(h^'^(r) • ey) — (Sa(r) • ey)(h^^(r) • ex) is instant angular frequency of 
the spin rotation. At is the time integration step, r is the reduced time related to real time t as 
r = fio'jMgt, where fiQ is permeability of free space, 7 is the gyromagnetic ratio and Mg is the 
saturation magnetization. The unitarity of spin rotation matrix 
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implies that numerical scheme adapted from Eq.(2) conserves the spin vector size independently of 
jt+At ^^^^/-j^^/ quadrature. 

The time demanding task represents summation of dipolar fields. The substantial reduction of 
time allows the use of the hierarchical approach [20] with a detailed summation of near and rough 
summation of far field contributions. The specific geometry of array leads to the natural spin block 
definition 

L2_i 

S,^= E je {0,1,2,..., Ll-1}, (4) 

0=0 

where sum runs over the spins of index (3 belonging to dot index = j and dot position rf^^y In 
frame of two-level hierarchical approximation dipolar field is splitted to near h^^'^^ and far h^^'^ field 

K ' =-f 2^ r— ^ + 2^ e„/3S/3 + -r , (5) 
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depending on the unit vectors n^^^x . = rfr^x /Irf, x I, r'^r^^- = — rfr^x. The summation rule is 
controlled by the scale-separation cut-off rcut- 

3 The results of array reversal simulations 

In this section we present the simulations of the square shaped finite array including 40000 spins. 
To analyze the simulation outputs we applied: (i) the qualitative recognition of the configuration 
snapshots; (ii) the comparing of the complex array dynamics with the dynamics of much simpler 
systems; (iii) the clustering of intra-dot magnetic configurations using neural networks. The global 
information on the array ordering concerns the magnetization per site 

M(r) = ^ E Sa(r). (7) 

The simulations were performed for uniform zero Sq;|t-=o = — H/iJ = ex and non-zero bias Sq.|,-=o = 
cos(99o)6x + sin(99o)ey, 9^0 7^ initial conditions. In both cases the characteristic time of array is the 
time when H • M projection gains at first the bound 

H-M(Teh)=0. (8) 



3.1 The qualitative approach based on two-scale similarity 

To obtain the preliminary qualitative information about the switching mechanism we run the simu- 
lation of reversal dynamics in the field H/J = — 1 and coupling Z)/J = 0.1 for dot size Lq = 20 and 
bias (po = 0. Under these conditions we obtained estimate t^i = 8.516. 

At the beginning of reversal the external magnetic field nucleates reversal domains at the corners 
of dots. The recognition of simulated systems shows that remagnetization is more homogeneous and 
easier when dots are localized at the corners of array, whereas nonuniformities occur due to mirror 
symmetry of effective fields typical for dots standing in positions closer to axial zone of array. This 
property is illustrated by snapshot taken at characteristic time as is shown in Fig. 1. The typical for 
nonuniform remagnetization is that reverse domains, which gain external field polarity, grow from the 
dot corners and pass toward the dot centers. At the intermediate stages domain walls - nonlinear waves 
[5] between reverse and initial polarity domains become narrow. In the context with nanoparticles 



and our small-scale simulations we prefer to use the term nonlinear wave instead of domain wall since 
nonlinear wave better expresses nonlocal character of arising spin structures. At the final stage of the 
inhomogcncous reversal the sum of external and dipolar field can, if sufficiently strong, homogenize 
intra-dot ordering by annihilating of wave/anti-wave pairs (WAP). Similar mechanism was studied in 
detail for linear spin chain [21]. As we see from the scheme depicted in Fig. 2 the reversal mechanism 
in dots is position and time dependent for finite size arrays. The formation of WAP is an example of 
inhomogcncous rcmagnctization. The process is characteristic for later times and dots placed nearly 
the axial array zone as is depicted in Fig. 3. The parallelism between spin chain and array behavior 
inspired the monitoring of spin projections summarized in Fig. 4. The figure shows Say components 
taken from several superlatticc site lines crossing array in external field direction. The time sequence of 
these projections reveals that each inter-dot interface acts as a phase-shifter of WAP formed in different 
dots. The consequence of inter-dot coupling is the large-scale WAP with shape resembling its small- 
scale intra-dot counterpart (see Fig. 5). It is rather surprising that for its build up magnctostatic forces 
suffice. The spatio-temporal two-scale similarity of magnetic ordering illustrated in Fig. 6 originates 
from spatial distribution of dipolar fields of square-shaped dots embedded into square-shaped array. 

3.2 The quantitative analysis of switching, the effect of bieis 

3.2.1 The dot-size dependences 

We continue with the comparing of Tch(-t'o) dependences for patterned structure as well as for the 
isolated dots with zero and non-zero bias. The varying of Lq is in some sense equivalent to varying 
of dipolar and exchange coupling contributions. The results for patterning Lq = 5,10,20,25,40,50 
for D/J = 1 are depicted in Fig. 7, where the characteristic times of the interacting and isolated 
dots are also compared. From the simulations it follows that inhomogcncitics grow as Lq increases. 
The reverse phase nucleates preferentially at dot surfaces and reduces Tch in larger dots. If inter-dot 
interactions are considered rch(i^o) dependences look smoother than dependences obtained for isolated 
dots. The interactions affect spatial dispersion of Wq,, which means that some dots arc accelerated 
and some decelerated during reversal. The simulations of isolated dots show that differences between 
zero and non-zero bias are larger for smaller Lq. The decrease of Tch due to uniform bias (^o > is a 
manifestation of the Stoner-Wohlfarth coherent regime. 

3.2.2 The interplay of D, J, H, ipQ parameters, minimal four-spin model 

In this section minimal four-spin model and complete array model of magnetization reversal are 
discussed and compared. The minimal model consists of four spins arranged to 2 x 2 square segment 
of unit size. In the segment four nearest pairs arc coupled by the exchange interaction and all of the 
pairs by dipole-dipole interaction. Despite of substantial geometric truncation, all of the interactions 
considered for array take place in four-spin model. The simulation of the minimal model system 
provides r^j^ characteristics depicted in Fig. 8. They consist of two qualitatively different regimes 
separated by bifurcation point D = manifested by t^i^D) cusp. At several modes compete, the 
interaction and external fields are nearly balanced and frustrated spins need quite long time to select 
the winning mode. For D < the reversal is very susceptible to choice of initial bias. It is easy to 
check that even small ipo > accelerates reversal and stabilizes coherent or quasi-coherent rotation in 
coincidence with the Stoner-Wohlfarth scenario. Further simplification of four-spin dynamics which 
clarifies situation can be done assuming the total coherence of spins. The consequence is the elimination 
of three angular variables and analytic expression 



where £ = AttX/juqMs is undimcnsional parameter proportional to damping parameter A. The 
expression obtained in the limit D ^ is rather well known 



For D > the integral Eq.(9) leads to finite r^j^ even in the case (po = 0. It means that dipolar 
coupling removes logarithmic divergence oc ImpQ of t^-^ and substitutes in part the effect of tpo > 0. 
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In the regime D < D^. the bias (pQ = induces non-uniform flower mode belonging to characteristic 
time which is out of the scope of Eq. (9) . The four-spin dynamics attained by simulation for D > 
is even more complex. In that case the choice of bias losses its dynamical relevance. The simulations 
show that quite short time is needed to form the transient long-living flower mode (long with respect 
to r^jj) which alters to vortex mode at final stages. 

With the four-spin system behavior in mind the attention is turned now to the complex array 
dynamics. As we see from Fig. 9 Tc\i{D / J) dependence does not contain sharp cusp. Instead of it, 
the slope variation of Tch{D/ J) arises (labeled as Dc in analogy to D^). For Lq = 20 and H/J = 1 
the position of Dc can be estimated as Dc — 5 J. The configuration snapshots show that near the 
array consists of spin configuration loops of diameters approaching single dot size or size of several 
dots. For D > Dc the array reversal generates non-collinear antiferromagnetic phases. The common 
feature of oversimplified four-spin model and array simulation is that finite t^-^ and Teh are attained 
even for tpo = due to non-collinearity of external and dipolar fields. From dependences depicted in 
Figs. 7 and 9 it follows that geometric variations in Lq can produce less remarkable changes than D/J 
parameters. 

For completeness the external field dependences of characteristic time should be mentioned. From 
Eq.(9) one can obtain for D/H ^ 1 the asymptotic expansion of characteristic time t^^i — {J/i^^)) 
Y,'^^Q{-lY{D/HYcn{^Q), Cni^o) = 3"2-W2 jV2 pQgn^2(^)(sin((/?))-'^-i suitable for proposition 
of reasonable fitting function covering the array data. The array simulations in Fig. 10 shows that Teh 
coincides qualitatively with asymptotic predictions of four-spin system. 



3.3 The intra-dot complex structures in representation of ART network 

The artificial neural network models [10] inspired by the physiology are able to mimic action of neurons 
and synaptic connections of brain. The salient feature of neural systems is the associative recognition 
of the complex structures. This ability is attained by training of synaptic connections - neurons. 
In this paper we deal with networks based on the adaptive resonance theory (ART) developed by 
Carpenter and Grossberg [22] originally as a model explaining adaptive phenomena in visual systems. 
The family of ART algorithms belongs to unsupervised self-organized algorithms motivated by need to 
construct network sufficiently adaptive to novel inputs. The adaptivity means that network produces 
new neurons if diversity of inputs overcomes the predefined threshold represented by so called vigilance 
parameter. On the other hand redundancy of information leads to the annihilation of neurons. For 
summary of ART network applications see e.g. [23] 



3.3.1 The encoding of intra-dot magnetic structure 

The ideas of clustering and adaptivity of the ART networks led us to their application to magnetic 
intra-dot configurations generated during reversal. At first, an appropriate format for encoding of 
magnetic configurations is suggested. Similarly to formulation [11] the configurations are encoded by 
2 A^c dimensional tuples 

mi = [ma,mi2, . . . ,miArJ , (11) 
where 

^in = ^J2^l^ i = l,2,...,L^, n = l,2,...,Nc (12) 

is the local average taken from Ns spins (see Fig. 11), where sum runs over site / taken from nth 
segment of ith dot labeled as O^^. Each dot is subdivided to Ac = Lq/Nq, intra-dot square segments. 
The system of intra-dot configurations is transformed into tuples incoming to ART, which 
compresses them to < L| neurons 

Wj = [wji,Wj2, WjiVc] , j G {1, 2, . . . N^} (13) 

representing the basal types of intra-dot ordering. The self-organization mechanism of neurons dis- 
cussed in the next is controlled by the coincidence between input rhi and output Wj. As a proper 
measure of coincidence is chosen Euclidean distance 



\W-i — m,; 



\ '^i^jn - mjn) • (Wj„ - mj„)T , (14) 
\ n=l 



where T labels transposition. 



3.3.2 The ART clustering algorithm 

The training algorithm is described in the following points 

1. Initialization setting N„ = 1, wi = rhi for random dot selection between i £ {0, 1, . . . , — 1}. 

2. Loop through the set of intra-dot magnetic configurations. For each randomly selected dot i 
follows: 

2.1 Presentation of rhi to ART network. The training is realized for k = 1,2,... iteration steps. 

2.2 Computing of the actual index j*(k, i) of the fired neuron ii'j*(fc,i)(fe) trained by rhi pattern 
according to the competitive rule 

j*{k,i)= arg min \\wj{k)-rhi\\ , . 

j=l,2,...,iVw • ilOj 

2.3 Comparing — rhi\\ to the vigilance parameter p. 

2.3.1 Update of the neuron weights via Hebbian training rule [10] 

Wj*{k,'i){k + 1) = Wj*(k,i){^) + V{k) (rhi - Wj.(^k,i){k)) (16) 

is applied if \\'Wj*(^k,i) ~irii\\ < P- The update shifts winning vector toward nearby input 
with rate proportional to plasticity parameter r]{k) = ?7o exp (— fe/rtrain) relaxing with 
training time rtrain- 

2.3.2 Creation of new neuron ^ + 1, Wj^^ <— rhi if ~ "^ill > P- 

2.4 ^nm/iiZation of neuron pair ['Wz*ik),Wz*{k)], ordered as < with indices , Z2 selected 
according to minimum distance 

[4^4] = arg min - WziW n?^ 

if \\wz^ — Wz* II < p • The product of annihilation [2;^, -Z2] ~^ -^w — 1 is the neuron 

Wz* determined by the midpoint rule Wz* ^{wz* + Wz*)- After this update of Wz* the 
neurons having j > Z2 are rewritten: Wj-i Wj. 

3. Stop criterion is represented by the inequality 

1 ^^-^ 

72 Yl II*j*(m(^ + 1) - '^j*{k,i)i^)\\ < £ ' (18) 

where e is small parameter. If the above inequality is not fulfilled, the algorithm follows from 
the step 2 with k incremented by 1. 

The stop criterion indicates when the network attains a fixed point of Wj, j = 1,2, . . . including 
the cluster of typical representative intra-dot configurations. 



3.3.3 The ART applied to different stages of reversal 

The ART network has been applied independently to array configurations attained at selected stages 
of reversal. The results of self-organization are presented in Figs. 13-15. The ART self-organization 
was carried out for tuned parameters rjo = 0.1, Ttj-ain 

= 20, £ = 10 ^, p = 1.15. The most important 
seems to be optimized choice of an appropriate vigilance connecting p and iV^. It stem from calibration 
in Fig. 12. During it p was decreased slowly from the initial p = 1.5 stabilizing single neuron. This 
decrease of p was stopped for A/^ = = 100 competing neurons. 

The ART view point on reversal leads to intra-dot taxonomy depicted in Fig. 13. In this figure the 
most similar neurons [in the sense of generalized distance Eq.(14)] at different times are joined by the 
arrows. The ART model clearly distinguishes between branches of coherently rotating monodomain 
dots and dots with more complex order. The diversity of configurations occurring near the internal 



energy maxima is reflected by the cardinality of neuron population. The ART pruning at late time 
stages is associated with small intra-dot diversity within the array. How choice of p affects the 
information content of network is clearly demonstrated in Fig. 14. The supplementary problem solved 
by ART represents classification of interfacial inter-dot configurations. In that case the configuration 
tuples m,, Wj have been combined from nij^ segments [see Eq.(12)] of adjacent quadrants of four 
neighboring dots. The results depicted in Fig. 15 indicate preferential formation of inter-dot vortices 
and strong tendency to inter-dot pinning also for conditions when formation of intra-dot vortices is 
disadvantageous. 



4 Conclusions 

The magnetization reversal has been simulated for superlattice spin model of the ultra-dense magnetic 
dot array. Several aspects of reversal dynamics for different exchange and magnetostatic couplings, 
external fields, geometric parameters and initial conditions have been studied. It was found that 
the principal physical consequence of geometric modulation is the formation of two-scale WAP. The 
ART viewpoint seems to be fruitful in thinking how dot switches from one configuration to another 
within the large time intervals, and how to identify dominant channels of intra-dot evolution. The 
present experience creates believe that advantages of ART paradigm should be exploited namely in 
3-D lattices, where visual classification meets natural bounds. 
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A Appendix 

A.l One dimensional enumeration scheme of superlattice sites 

In the Appendix site enumeration scheme mentioned in section 2 is described. In this scheme the 
superlattice relative position vector Fq,^ pointing from site a to site (3 is equal to 

ra/3 = [M{P,L)- M{a, L)]e^ + L) - X(a, L) ] ey , (19) 

where 

M{a,h) = a mod 6, I{a,b) = [a/b].^^ . (20) 

Here Ai is standard modulo function and I is the integer part of its argument {L is an integer multiple 
of Lq). The second part of Eq.(l) includes the exchange coupling Je^/j generated by dependence 

ea/3 = ( 1 - S<t>{l3)-l,<l>{a)-Lo) ( 1 " '^^(/3)-l,^(a)-Lo ) " (21) 

It gives rise to geometry of dense square patterning of the period Lq in perpendicular ex, By directions. 
In Eq.(21) 6 is the Kronecker symbol and 

</.(«) = M{M{a, L), Lo) , C(a) = M{I{a, L), Lq) (22) 

are dependences illustrated in Fig. 16. The position of some dot including site a can be also written 
in terms of M and I functions: 

rf(„) =I{M{a,L),Lo)e^ + I{I{a,L),Lo)ey. (23) 

Here 

i{a)=L^I{I{a,L),LQ)+I{M{a,L),Lo) (24) 
is one-dimensional map associating site index a with dot index i{a) G {0, 1, 2, . . . , — 1}. 



A. 2 Integration scheme 

To derive explicit integration scheme with variable time steps KiAr and k,2^t it is supposed that 
integral of spin angular frequency uJq, (r) is approximated by 



t+At 

a;„(T')dr' = AraJa(r) , (25) 



where 

^aij) = «0 ^aij) + 0-1 ^a{T ~ Kl Ar) + 02 (^^(t — Ki Ar — K2 At) (26) 

is r— local Ar averaged Uair') expressed through coefficients 

2 + Ck? +3k2 + 6/Vl(l + K2) 

6«;i(ki + K2) 
2 + 3K1 + 3k2 

ai = 7 , (28) 

bKlK2 

2 + 3K1 

= 7—7 T ^ (29) 

0K2{Ki + K2) 

derived by equating to zero the series of difference [lhs.Eq.(25) — rhs.Eq.(25)] considered up to the 

order ©(Ar^). The reason for use of non-constant prefactors ki,2 ~ 1 is the opportunity to change 
adaptively integration step via the instantaneous spin rotation extreme 

max uJair-KiAr). 
a=0,...,L2-l (^0) 



In the uniform case ki = K2 = 1, the quadrature is compatible with Adams-Bashworth 3th order 
integration formula. The initial two steps estimating LOaiT — Ar) and u>{t — 2 At) can be performed 
using Euler scheme obtained for ao = 1, ai = 02 = 0. 
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B Figure Captions 




Figure 1: The array snapshot including configuration with the resolution given by A'^s = 4 = 2x2. 
Obtained for parameters D/J = 0.1, Lq = 20 at characteristic time Teh = 8.516. The differences 
between coherent rotation at array corners and flower mode near to the array center are striking. The 
dashed lines inside magnified segments a), b) belong to bounds separating dots. Along these lines the 
exchange couplings are broken. 





time 



Figure 2: Two schemes showing principal transients of nucleation, WAP formation and annihilation. 
The intra-dot open (straight) WAP (case a) occurring at earlier transients of remagnetization in 
eccentric dots. The center-symmetric WAP closed mode (case b) typical for late transients and dots 
positioned near the array center. 




Figure 3: The array snapshot and its magnified details a), b) obtained for late stage (r = 11) of 
reversal (the dot labels and model parameters are identical as in the Fig. 1) nearly the time when the 
internal energy (free from Zeeman term) reaches its maximum. 
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Figure 4: The reversal in the representation of one dimensional spin projections ^i{t) ■ Cy. The left 
hand side sequence of figures constructed for 22'^^^ stright line (solid line) and for 37*^^ line (dashed 
line) in the external field direction ex- The right hand side sequence belongs to 97*^ (solid line) and 
102°'^ line (dashed line) of the considered system D/J = 0.1, Lq = 20. 




Figure 5: The snapshot illustrating similarity of nonlinear WAP at intra-dot (Lq = 20) and inter-dot 
(L = 200) spatial scales (at different times). The rescaling is performed for dot indexed by i = 55 taken 
from 102"'^ line of array [spin index changes within a G (20500, 20519)] at time r = 13.5 (depicted 
by solid line) and for 28*^ spin line of array [where a G (5600,5799)] at r = 9.7 (dashed line). The 
intcr-dot and intra-dot configurations look similar when the spin coordinates are transformed using 
multiplicative factor Lq/L. 



Figure 6: The spatio-temporal similarity of the magnetic configurations occurring at intra-dot and 
array scales. On the left hand side we plot overall moments of dots at different times r = 4.5 (upper 
array flower) and r = 8.5 (bottom array configuration). The detailed intra-dot configurations shown at 
the right hand side column corresponding to dot moments marked by the bold vector and fixed at the 
later times r = 8.5 (upper intra-dot flower) r = 12.5 (bottom three-domain intra-dot configuration) 
for D/J = 0.1, Lo = 20, L = 200. 




Figure 7: The Tch(-^^o) dependence calculated for D/J = 1, H/J = —ex- The insert compares 
characteristic time for ipo > and ipo = configurations in isolated dot. The most remarkable are 
differences for small Lq. The lines are guidelines to eye. Only a weak dependence on hierarchical 
cut-off is observed. The test of T^h accuracy for a given size of array is performed for Lq = 5 and 
^cut £ {V^Lq,2\/^Lq,3\/2Lq}. The calculations confirm that changes in Teh does not exceed bound 
0.004. 
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Figure 8: The Tch{D/ J) dependences obtained for different remagnetization (flower and vortex) modes 
of four-spin system simulated for ipQ > (left hand side) and for (pQ = (right hand side). The cusp 
corresponds to four-spin bifurcation point D^. 




Figure 9: The plot of Tch{D/J) dependence obtained for modulation Lq = 20. The variation of 
TchiD/J) slope near Dc associated with change of reversal mode. 




Figure 10: The plot of Tch{H/J) dependence obtained for Lq = 20 and D/J = 0.1. Dashed hne 
represents fit Teh = a{J/H) - bD {J/H"^) with parameters a = 8.302 ± 0.122 and b = 1.506 ± 0.131. 
The choice of fitting function stems from the asymptotic iJ-dependence derived for four-spin model 
in sect. 3. 2. 2. 
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Figure 11: The schematic example of intra-dot segmentation. Figure showing sixth square segment 
labeled Dje belonging to ith dot including iVs = 49 spins. 
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Figure 12: The pruning of ART neurons. The number of neurons decreases with increasing vigi- 
lance parameter. The dependences obtained for D/J = 0.1 and Lq = 20. The intra-dot configuration 
is locally averaged over A'^s = 4 spins [see Eq.(12)]. Case a) corresponds to configuration near to the 
switching time (r = 8.5). Case b) belongs to time when the internal energy (free from Zeeman term) 
is near its maximum (r = 11). 
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Figure 13: The ART network representation of the main intra-dot transients of array reversal for 
D/J = 0.1, Lq = 20. The configurations corresponding to different times obtained for As = 4 
segment average [see Eq.(12)]. The arrows connect closest transient configurations (in the sense of 
the magnetic Euclidean distance). The extremal number of neurons occurring for r = 9.1 roughly 
corresponds to the maximum of internal energy. 
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Figure 14: Figure showing ART trained with relatively small vigilance p = 0.8. The network properly 
uncovers mechanism of WAP annihilation and nonlinear wave motion. 
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Figure 15: The ART revealing interfacial four-dot modes of magnetization reversal: inter-dot nonlinear 
wave pinning and formation of inter-dot vortices. Parameters: D/ J = 0.1, p = 1.22, A'^g = 4, Lq = 10. 
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Figure 16: The illustration of one-dimensional site enumeration scheme (Lq = 3, L = 9). The 
modulation of connectivity matrix e^/j with period Lq = 3 in ex and Cy perpendicular directions 
realized by functions 0(a), ^(a) from Eq.(21). 



